Non-equilibrium fluctuations and mechanochemical couplings of a molecular motor 
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We investigate theoretically the violations of Einstein and Onsager relations, and the efficiency 
for a single processive motor operating far from equilibrium using an extension of the two-state 
model introduced by Kafri et al. [Biophys. J. 86, 3373 (2004)]. With the aid of the Fluctuation 
Theorem, we analyze the general features of these violations and this efficiency and link them to 
mechanochemical couplings of motors. In particular, an analysis of the experimental data of kinesin 
using our framework leads to interesting predictions that may serve as a guide for future experiments. 
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Motor proteins are nano-machines that convert chem- 
ical energy into mechanical work and motion Im- 
portant examples include kinesin, myosin, and RNA 
polymerase. Despite a number of theoretical models 
[3, 13, 0, S @| , understanding the mechanochemical trans- 
duction mechanisms behind these motors remains a sig- 
nificant challenge Recent advances in experimental 
techniques to probe the fluctuations of single motors 
provide ways to gain insight into their kinetic pathways 
[10|. However, a general description for fluctuations of 
systems driven out of equilibrium, and in particular of 
motors, is still lacking. Recently, the Fluctuation Theo- 
rem (FT) [111, [13, [III has emerged as a promising frame- 
work to characterize fluctuations in far-from-cquilibrium 
regimes where Einstein and Onsager relations no longer 
hold In a nutshell, FT states that the probabil- 

ity distribution for the entropy production rate obeys a 
symmetry relation, and it has been verified in a number 
of beautiful experiments on biopolymers and colloidal 
systems In this Letter, we demonstrate that FT 

provides a natural framework in which thermodynamic 
constraints can be imposed on the operation of nano- 
machines far from equilibrium. 

Specifically, we study a generalization of the two-state 
model of motors introduced in Ref. [5|]. Although simi- 
lar models have been investigated with known exact re- 
sults we reformulate the model to include an im- 
portant variable, namely the number of ATP consumed, 
and construct a thermodynamic framework. Our frame- 
work allows us to characterize the ATP consumption rate 
of a motor, its run length, and its thermodynamic effi- 
ciency. Additionally, we show that our model obeys FT 
p^ . While there have been a few recent studies prov- 
ing FT for motors [l^, [l^, [l3| , we further investigate the 
physical implications of FT here. In particular, we quan- 
tify the violations of Einstein and Onsager relations, re- 
spectively, by four temperature-like parameters, , and 
by the difference of the mechanochemical coupling coeffi- 
cients, AA, and we explore the behaviors of and AA, 
as well as the motor efficiency, as functions of generalized 
forces with the aid of FT. Our main results are (i) one of 



the Einstein relations holds near stalling, (ii) the degree 
by which the Onsager symmetry is broken (AA ^ 0) is 
largely determined by the underlying asymmetry of the 
substrate, (iii) only two "effective" temperatures char- 
acterize the fluctuations of tightly coupled motors, and 
(iv) kinesin's maximum efficiency and its maximum vio- 
lation of Onsager symmetry occur roughly at the same 
energy scale, corresponding to that of an ATP hydrolysis 
(- 20 ksT). 

As a result of conformational changes powered by hy- 
drolysis of ATP, a linear processive motor, like kinesin, 
moves along a one-dimensional substrate (microtubules). 
Its state may be characterized by two variables: its po- 
sition and the number of ATP consumed. To model its 
dynamics, we consider a linear discrete lattice, where the 
motor "hops" from one site to neighboring sites, either 
consuming or producing ATP (see Fig. [T|) . The position 
is denoted hy x = nd, where 2d « 8 nm is the step size 
for kinesin. The even sites (denoted by a) are the low- 
energy state of the motor, whereas the odd sites (denoted 
by b) are its high-energy state; their energy difference is 
AE = ksTe, where is the Boltzmann constant and T 
is the temperature. 

Because of the periodicity of the filament, all the even 
(a) sites and all the odd (6) sites are equivalent. The dy- 
namics is governed by a master equation for the proba- 
bility, Pi{n, y, t), that the motor, at time i, has consumed 
y units of ATP and is at site i {— a,b) with position n: 

dtPiin,y,t) = - (tJi -I- Toi) Pi{n,y,t) + 

E 

;=-i,o,i 



[ t7j Pj(n +l,y-l,t) + ^] Pjin -\,y-l,t)' 



with i ^ i, where oj j — and lo j = ujj. 

Denoted by tjj and To ^ are the transition probability 
per unit time for the motor, with l{— —1,0,1) ATP 
molecules consumed, to jump from site j to a neighboring 
site to the left or to the right, respectively. 

The transition rates can be constructed by consider- 
ing the kinetics of the transitions between the two states 
Mq and Mft of the motor f3l|. We assume two different 



FIG. 1: A schematic of the rates of the two- state model for 
a molecular motor moving on a linear lattice with y number 
of ATP consumed. The even and odd sites are denoted by a 
and 6, respectively. In the case of two headed kinesin, site a 
represents a state where both heads are bound to the filament, 
whereas site h represents a state with only one head bound. 

chemical pathways: (a) + ATP ^ Mb + ADP + P, 
{(3) Ma ^ Mfe. The a-pathway represents the transition 
of the motor accompanied by ATP hydrolysis and the [3- 
pathway represents the transition driven by thermal acti- 
vation. It is straightforward to generalize the model with 
more chemical pathways, but here we focus only on these 
two, for which uTb^ = uJ^~^ = uj^~^ — ujI^ = 0. Following 
Ref. Q , the transition rates in the presence of an external 
force Fe are changed according to tJ'-{Fe) — tj -'(0) e^^* 

and a;/(Fe) = w/(0)e+^^/, where / = Fed/ikeT) and 
of are the load distribution factors Q. These load dis- 
tribution factors take into account that the external force 
may not distribute uniformly among different transitions. 
After one period, the work done by F^ on the motor is 
-Fe2d, implying that 0+ + e^ +9- +6+ = 2 [|. Thus, 
we may write the non-zero rates as: 

uJl^^ — a e^^b ^ , LOiP = LD e~^b ^ , 

where a and a', and to and u' are the bare rates for 
the two distinct transitions for the pathways and A/I = 
kBTAfi is the chemical potential difference f3|. The un- 
derlying asymmetry of the substrate dictates that a ^ a' 
and w 7^ cj' as required for directional motion [5]. 

We find that the rates in Eq. ([T]) satisfy four generalized 
detailed balance conditions: 

^b-'Ph"" = U-Jp-'ie+^'^+'^'>f-^^^\ (2) 
Ul-ip^'i = ^JP^^e-('i+<'^^f-^^\ (3) 

for / = 0,1, where P^'i = 1/(1 + e"') and P°'^ = 
6^*^/(1 -|- e"*^) are the equilibrium probabilities corre- 
sponding to / = and A/i = 0. We note that these 
relations, Eqs. ([2]) and ([3]), while valid arbitrary far from 
equilibrium, still refer to the equilibrium state via the 
probabilities P^. We show below that these relations 
lead to a FT |l3|- Introducing the generating functions: 
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FIG. 2: Four modes of operation of a molecular motor, as 
delimited by t) = and r = Ql- The lines are generated 
with parameters that we have extracted from fitting the data 
for kinesin in Ref. 0| to our model. (The best-fit values for 
the parameters are listed at the end of the text.) In Region 
A, where rAfi > and fv < 0, the motor uses chemical 
energy of ATP to perform mechanical work. In Region B, 
where rA/i < and fv > 0, the motor produces ATP from 
mechanical work. In Region C, where rA^ > and fv < 0, 
the motor uses ADP to perform mechanical work. In Region 
D, where rA/i < and fit > 0, the motor produces ADP 
from mechanical work. 

F^{zl,Z2,t) = EyEne~"''"~"'''P^i'^^y^^)^ ^hosc time 
evolution is governed by: dtFi = Mij Fj, where [^i , Z2] 
is a 2 X 2 matrix that can be obtained from the master 
equation above, we find ( e^^i"^^^?' ) = ^ . Fi{zi,Z2, t) ^ 
exp(i?i), for t — > 00, where d = i?[zi,Z2] is the largest 
eigenvalue of M. Using Eqs. ^ and ([3|), it can be shown 
that A4 and are related by a similarity transforma- 
tion: M^[f - zi,A^i- Z2] = QM[zi,z2] Q"\ where 
is the adjoint of A4 and Q is a diagonal matrix. This 
similarity relation implies that 

d[z,,Z2]^^[f-Z,,AfX-Z2], (4) 

which is one form of FT. 

Now, we proceed to discuss the physical consequences 
of FT. The eigenvalue, d, contains all the steady-state 
properties of the motor. In particular, the average 
(normalized) velocity, v = w/d, and the average ATP 
consumption rate, r, are, by definition, given hy v = 
~d,,i}[0,0] and r = -d^^'d[0,0], respectively [Hi. From 
the conditions: v — and r = 0, we can construct a full 
operation diagram of a motor, as shown in Fig. [2] for the 
case of kinesin. The curves v = and r = define implic- 
itly / = /st(A/i) (the stalling force) and Afi = Afi^tif), 
respectively. It is interesting to note that the large asym- 
metry between regions A and C in Fig. [2] reflects the fact 
that kinesin is a unidirectional motor. 

The response and fluctuations of a motor are quan- 
tified, respectively, by a response matrix Ay and by a 
diffusion matrix: 2Dij = 9zi9zj i?[0, 0]. The physi- 
cal meanings of are: An = dv/df is the mobility, 
A22 = dr/dAfi is the chemical admittance, and more 
importantly, A12 = dv/dAfj, and A21 = dr/df are the 
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Onsager coefEcients that quantify the mechanochemical 
couphngs of the motor. Differentiating Eq. we can 
write: 

Near equihbrium, where / and A/i are small, a Taylor 
expansion of Eq. (O leads to w = A^^ / + A/i and 
r = X21 f + Aa*, with A°- = d^^d,^/[0, 0] = Aj, which 
are the Einstein relations, and A°2 = dz2dzi'&[0,0] — 
dzidz2^[0,0] = A21, which is the Onsager relation. Thus, 
FT captures the response and fluctuations near equilib- 
rium (l^ . [l7| . 

Away from equilibrium, we expect that Onsager and 
Einstein relations are no longer valid. To quantify 
their violations, we introduce AA = A12 — A21 and 
four "temperature" -like quantities, Tij = Dij/Xij. Of 
course, these effective temperatures are not thermody- 
namic temperatures: they are merely one of the ways 
to quantify deviations of Einstein relations. Via FT, we 
obtain the following general characterizations for these 
quantities. First, at sufficiently small driving, we find 
2AA « (5a^^9ii - df{)i2)f + (9a,.i?i2 - 5/i?22)A/x, where 
•dij = dzidzj /2, Afi/2]. In particular, for / ^ 1, 
AA cx Afi. Thus, active processes in the mechanochem- 
ical transduction mechanism breaks Onsager symmetry 
and these processes can be studied via AA. Secondly, 
along v{f, A/i) — 0, we find that Eq. ^ has a special re- 
latiom d[zi,0] = d[6f - zi,0], where Sf = f - U{A^i). 
Therefore, one of the Einstein relations. An = -Dii, holds 
near stalling, since FT implies that 2v = 0,0] 6 f for 
small Sf. Note that this particular Einstein relation also 
holds for ratchet models under similar consitions [19|] . By 
the same token, near r = 0, FT implies that an Einstein 
relation holds for y, i.e. D22 — X22- 

For our two-state model, we can fully investigate the 
behaviors of AA and Tij. Let us focus on region A of 
Fig. [5] and — / ^ 1, so that Ay and Tij depend only on 
A/i. In Fig. [3^, we display AA and the three distinct 
Tij (see below) as a function of Afi. We observe that for 
small A/i, AA rises linearly with A/x, in agreement with 
the FT prediction, and that for larger A/i, AA exhibits 
a maximum. Moreover, for large A/i, we find that AA 
approaches to a constant value. The latter observation 
can be understood from a simple argument. When A/i ^ 
1, the transitions between the states a and b of the motor 
are limited by the /3-pathways. Therefore, we can write 
r ~ u!e~^i> ^ + ui' e^i> ^ , which implies that for small /, 
AA ^ uj — Lu' 9'^ since A12 « for large A/t. Thus, the 
underlying asymmetry of the substrate determines the 
degree by which the Onsager symmetry is broken. 

In Region A and — / ^ 1, the also exhibit interest- 
ing behaviors. First, we note that the run length £ - the 
distance moved per ATP hydrolyzed - is independent of 
A/i: e = v/r = 2d {a lu' - a'to) /[{a a') (w + w')] < 2d. 
With the help of FT, we find that T12 = T22 for any 
A/i. Therefore, there are only three "effective" tempera- 
tures instead of four, as one might naturally suppose. As 
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FIG. 3: (a) Plots of Tn (dot-dashed), T21 (dotted), T22 
(dashed), and AA (solid) vs. A/t in Region A of Fig. [2] with 
small /. Tij characterize the fluctuation-response ratios (see 
text), and AA quantifies the breaking of Onsager symmetry, 
(b) Local maximum of the efficiency 77m vs. A/i. Note that 
77m is substantially larger than rj^ (the dotted line). Note also 
that the absolute maximum, which occurs at about A/i ^ 15, 
roughly corresponds also to the maximum of A A. Inset: Effi- 
ciency vs. normalized force for A/t — 15. The parameters used 
in both (a) and (b) are the same as those used to generate 
Fig. [3 



shown in Fig. [3^, we observe that all distinct Tij start 
off at Tij = 1 near equilibrium, as expected, and for 
large A/i, T22 ^ e^^ diverges exponentially whereas T21 
and Tn approach to finite values. Secondly, for tightly 
coupled motors, £/{2d) ^ 1, we find that Tn is nearly 
identical to T21 (see Fig. [3Ji). Therefore, in this case, 
only two "effective" temperatures characterize motors' 
fluctuations. 

In addition, our framework allows us to investigate the 
thermodynamic efficiency, an important quantity that 
also characterizes the working of a motor [23] . In re- 
gion A, it is defined as the ratio of the work performed 
to the chemical energy input: 77 = — /^/(rA/t) By 
definition, 77 vanishes at / = and at the stalling force 
fst- Therefore, it has a local maximum 7/m(A/i) for some 
/m(A/t) between fst < fm < (see Fig. inset). Near 
equilibrium, 77m (A/t) has a constant value, ttJ^, along a 
straight line /in(A/t) cx A/i inside region A [2!]. Far from 
equilibrium, we find that r/m has an absolute maximum 
at some A/i > 1, and rj^ is substantially larger than rj'^'} 
as shown in Fig. [5Jd. Hence, a motor achieves a higher 
efficiency in the far-from-equilibrium regimes Q . 

Finally, to discuss the relevance of our framework to 
kinesin, we carried out a global fit of kinesin velocity data 
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FIG. 4: Kinesin velocity vs. ATP concentration under an ex- 
ternal force. The solid curves are the fits of our model to 
data from Ref. [^. From the top down, the plots are for 
Fs = —1.05,-3.59, and — 5.63pN, respectively. Inset: Ki- 
nesin velocity vs. force under a fixed ATP concentration. The 
solid curves are fits to the data of Ref. From the top down, 
the plots are for [ATP] = 2 mM and 5 fiM, respectively. 



of Ref. 8] to our model at different external forces and 
two curves of force vs. velocity at different ATP concen- 
trations (see Fig. 2]). Assuming that e'^'^ = feg [ATP], we 
obtain the best-fit values for the parameters: e = 10.81, 
ko = 1.4 • 10VM"\ a = 0.57 s-\ a' = 1.3 • lO^S"!, 
uj = 3.5 s-i, w' = 108.15 s-\ 9+ = 0.25, 6*" = 1.83, 
6+ = 0.08, and 6' = -0.16. These values are reason- 
able within the accepted biophysical picture of kinesin 
[l[. First, e and fcj^^ represent the typical binding en- 
ergy 10 ksT) of kinesin with microtubules and the 



ATP concentration at equilibrium (~ 10~^ fiM), respec- 
tively. Secondly, 9~ = 1.83 indicates that the back-steps 
(transitions a ^ &) of kinesin contain most of the dis- 
placement sensitivity [l|. Moreover, our framework al- 
lows us to estimate a maximum stalling force of — 7pN, 
and more importantly, a run length of £ ~ 0.97(2d) 
and a global ATP consumption rate of r ~ llls~^, 
all in excellent agreement with known values [^. Us- 
ing the above parameters, we constructed the diagram 
of operation for kinesin (Fig. we made predictions 
about AA and Tij (Fig. and we obtained the effi- 
ciency for kinesin (Fig. [Sb). In particular, we find that 
Til 10 T, the maximum value of AA ~ 45pN"^s"\ 
and AA — lOpN^^s"^ at large A^. Under typical 
physiological conditions (A/I ~ 10 — 25 ksT), kinesin op- 
erates at an efficiency in the range of 40 — 60%, also in 
agreement with experiments []|. Lastly, we point out a 
remarkable feature: the absolute maximum of occurs 
approximately at a A/i at which AA is also a maximum, 
corresponding to an energy scale of 15 — 20 ksT (see Fig. 
ini). It is interesting to note that kinesins operate most 
efficiently in an energy scale corresponding to the energy 
available from ATP hydrolysis. 

In conclusion, FT links a set of physical quantities that 
reveal the mechanochemical couplings of a motor and our 
results support a growing consensus that FT provides a 
possible organizing principle for driven active systems. 
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